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1. Introduction 

Modem probes of material properties, such as the new inelastic neutron scattering 
facilities, are reaching such unprecedented sensitivity that they can measure the 
spectrum not only of a single quasiparticle excitation, but even two-particle excitations 
(e.g. [1]). These quasiparticles can collide, scatter, or form bound states just like 
elementary particles in free space. The spectrum of the multiparticle excitations is a 
crucial indicator of the underlying dynamics of the system. 

The experiments measure scattering cross- sect ions, which are proportional to the 
appropriate 'structure factor' for the system or material at hand [21 E]. It is therefore 
of particular interest to explore the critical behaviour of these structure factors in the 
vicinity of a quantum phase transition. In this paper, we present a review of this topic, 
comparing the theoretical predictions with some exact analytic results and numerical 
calculations for various models. We concentrate here on quantum spin models, but the 
major conclusions are applicable more generally. 

The theoretical behaviour of the total structure factor has been discussed since 
early days. More recently, people have begun to discuss the breakdown of the total 
structure factor into its component multiparticle contributions from one, two, .. etc. 
intermediate quasiparticles. Sachdev [3], for instance, discusses the behaviour of the 
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1-particle structure factor in his book on quantum phase transitions. In Section [2] of the 
paper, we draw together these theoretical discussions. 

In the remainder of the paper, we review the behaviour of the structure factors 
for some specific models. In Section [3] we look at the transverse Ising chain, which is 
exactly solvable, and hence yields some exact results for the 1-particle structure factors 
[5]. In SectionlH we review some numerical results obtained by series expansion methods 
for some other models, namely the trasverse Ising model in higher dimensions [5]. the 
alternating Heisenberg chain P, [7] , and the bilayer Heisenberg antiferromagnet [8] . 

Our main conclusions, in Section [5], concern the relationship between the 1-particle 
structure factor and the total structure factor. It is usually assumed that the 1-particle 
term dominates the total structure factor, and their scaling behaviour is the same; 
but this is not always strictly true. In the transverse Ising model and the dimerized 
alternating chain, for example, te 1-particle structure factor actually vanishes at the 
critical coupling, everywhere except at the critical wavevector. Only for the bilayer 
model does the 1-particle structure factor remain dominant at all wavevectors. This 
latter behaviour, however, is presumably more typical in generic quantum spin systems. 



2. Review of Theory 

Assuming magnetic scattering from atomic spins Sj localized on sites i of a Bravais 
crystal lattice, the neutron scattering cross section can be directly related to the 
dynamical structure factor [2] 

5"^(k, uj) = --Y. dt e^(-*-k(rj-rO)^«7( _ t) (1) 

where 

C-\r^-T,,t)=<SJ{t)Sj{0)> . (2) 

Here i,j label sites of the lattice, 0,7 label Cartesian components of the spin operator 
S, is the number of lattice sites, C°'^(r, t) is the spin-spin correlation function, and 
the angular bracket denotes the thermal expectation value at finite T or, at T = 0, the 
ground-state expectation value. The structure factor satisfies the condition of 'detailed 
balance' 

5"^(k,a;) = e^"5'^"(-k, -uj). (3) 

where (3 = l/ksT in the exponent is the usual Boltzmann coefficient. The time 
dependence of the spin operator is given as usual by 

Sf{t) = e'^'Sf{0)e~'"' (4) 

Integrating over energy gives the 'integrated' or 'static' structure factor 

/oo 1 
dtu ^°^(k, uj) = —Y, e*'^-('''-''j) < SfS^ >, (5) 

the spatial Fourier transform of the 2-spin correlator at equal times. 
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Integrating over momentum then yields a sum rule: 

1 , . roo 1 ^ ^ 

]v E y_ du s-\k, ^) = ^ E < s:s^ >, (6) 

involving the expectation value of two spin operators at the same point. 

If S°' and S"^ are Hermitian conjugates, which is usually the case of most interest, 
we can introduce a complete set of energy eigenstates |n > in equation ([T]) and integrate 
over time to get 

S-''{k,u) = Y.Sr{k,u), (7) 

n 

i.e. a sum over 'exclusive' structure factors or 'spectral weights' S^'^, where 

S:^{k,u;) = Lj26{uj-E^ + Eo)\Y: < ^n\S^\^o > e^'^^'f (T = 0) (8) 



or for T 7^ 



S:^ik,u;) = 2_J2siuj-E^ + E^)e-^^-\J2< ^Pn\S^\^m > e^'^'^f (9) 



m.n 



where En is the energy of the nth eigenstate, \ipo > is the ground state, and Z is the 
partition function 

Z = J2e-^^". (10) 

n 

If the system exhibits well-defined quasiparticle excitations, the intermediate states n 
can be classified into 1-particle, 2-particle or many-particle states, each state making a 
non-negative contribution, so that the total structure factor is real and positive semi- 
definite. 

Following Sachdev [¥], we may also define the corresponding generalized 
susceptibility x"'''(k, u;) by a Fourier transform in imaginary time [it r) 

'ik,un)= rrfr^C"^(ri,r)e-^('^-^'~-"^) (11) 

Jo : 



where lOn = 27rnT, n integer, is the Matsubara frequency arising from periodic boundary 
conditions across the strip of width f3 in imaginary time. Then x°''''(k, cu) for real 
frequencies is obtained by an analytic continuation iuJn ^ uj + i6, where 5 is a positive 
infinitesimal. The dynamic susceptibility measures the response of the magnetization S*" 
to an external field coupled linearly to S'^, oscillating with wavevector k and frequency 
LJ. One can show [2] that satisfies the Kramers-Kronig relation 

J-oo TT \l — LU 

where P indicates the principal part. 

If S"" and are Hermitian conjugates, then a fluctuation-dissipation theorem 
connects the structure factor S"'^ to the imaginary part of the dynamic susceptibility 
PS]: 

5-(k,u;) = -^^^-^/m{x-(k,o.)} (13) 
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If S"" and S'' are themselves Hermitian, one can show, using spectral analysis as for 
S"'^ above, that 

X'^^*(k,^) = x"^(-k,-^) (14) 

If both conditions are true, i.e. a = 7 and S°' is Hermitian, then the diagonal 
susceptibihty obeys 

X""(k,.;) = x""(-k,a;) (15) 

and 

X""(k,-c.)=x""*(k,a;) (16) 

Thus Im{x°"^} is an odd function of 00, while Re{x°"^} is an even function of uj. From 
( |T3|) . the dynamic structure factor then satisfies 

5""(k, -uj) = e-^^5""(k, u) (17) 
2.1. Critical Behaviour near a Quantum Phase Transition 

Now let us suppose that a quantum spin model undergoes a quantum phase transition 
as a function of some coupling A at temperature T = 0. The critical behaviour of the 
integrated structure factor can be obtained from a heuristic argument as follows. In the 
continuum approximation near the critical point, equation ([6]) for the static structure 
factor reduces to 



5"^(k) = I d^r e**^-"- < 5"(r)5^(0) >o 



where d is the number of spatial dimensions. 

The oscillating factor exp(2k ■ r) will kill off the contributions from large distances 
unless it is compensated by a corresponding oscillation exp(— zko ■ r) in the correlation 
function. Then we can write 



5'"'^(k) = j d'^r e^'i""^(r) (19) 

where q = k — ko, and g(r) is a smooth function. Scaling theory [9], [4] then tells us that 
in the vicinity of the critical point 

^(r)~r-(^+^-2+'')/(r/0 (20) 

where ^ is the correlation length, and z is the dynamic critical exponent. Thus when 
k = ko, the 'critical wavevector', we have 

5"^(ko) = J d^r r-("'+^~2+'?)/(r/^) ~ ^2-.-^ j ^dy y-{d+z-2+n)f(^y^ (2I) 

where y = r/^. As the coupling A Ac, corresponding to a quantum phase transition, 
we expect 

e~|Ac-A|-^ (22) 

and hence 

5"^(ko) ~ |Ae- Ar^'-"-'')^ (23) 
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where y' = g^y, so that at the critical coupling we expect S°''^{'k) to scale like g~(2-2;-»?) 
at small q. 

For the 1-particle exclusive structure factor, we may paraphrase Sachdev's argument 
[1] as follows. Assuming relativistic invariance of the effective field theory (i.e. z = 1), 
which applies to many though not all models, the dynamic susceptibility in the vicinity 
of a quasiparticle pole is expected to have the form 

X^'iK ^) = 21.2 ^A^^ ^•^2 + ■ ■ ■ (25) 

where e is a positive infinitesimal, c the quasiparticle velocity, A is the quasiparticle 
energy gap, and A"""' is the "quasiparticle residue" . Then the dynamic structure factor 
is 

5"^(k,.;) = -/m{x"^(k,^)} (26) 

TT 

Let 

E{\i) = x/c2k2 + A2 (27) 

then from (!25l) . (!26l) and (l27ll we can write the dynamic structure factor for the 1-particle 
state 

S?^{Ku^) = ^^S{u;-E{k)) (28) 

and hence the static structure factor 

5r;(k)= ^^^i7(k,^) = ^^ (29) 

where A°''^{'k) is the residue function, which in general may be a function of k. Note 
that S(k, u) at T = vanishes for u < 0, from equation (IT^ . 

From renormalization group theory [9], the scaling dimensions of these quantities 
are expected to be [5] dim[x"'''] = —2 + rj and dim[A'^'>'] = rj, or in other words we expect 
near the critical point 

A"^{ko)r^\K-\r, (30) 

E{ko) ~ |Ae-Ar, (31) 

and hence 

5r;(ko) ~ |A, - Ar(l-^)^ (32) 

just as for the total structure factor (recall here z = 1). In many cases, the 1-particle 
contribution will dominate the structure factor, but this is not always true, as we shall 
see. 
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These behaviours may be encapsulated in a scahng form. Assuming once again 
relativistic invariance of the effective field theory near the critical point {z = 1), so that 
the quasiparticle excitation energy is given by equation ( 1271) . and the energy gap 

A = E(ko)~|Ae-Ar, (33) 

then following Sachdev [1] the structure factor at low temperatures to one side of the 
transition is expected to take the form 

= (34) 

where is a universal scaling function and Zs is a normalization constant depending 
on the microscopic model. In the 'quantum critical' regime, A/T 0. 

At zero temperature, we may choose A as the reference variable rather than T, and 
write 

S(k..) = jL.i,(^,-) (35) 
or integrating over u, 

^(k) = ) (36) 

where 

^'sC{) = ^J^J^'^sC{,u;'). (37) 

If the energy gap is zero, as in the presence of Goldstone bosons, an energy scale can be 
constructed from the spin-stiffness ps or the Josephson correlation length - we refer 
to Sachdev [Ij for details. 

3. Comparison with exact Results 

3.1. Transverse Ising model in one space dimension 

The transverse Ising chain model is exactly solvable, and expressions for the energy 
spectrum, magnetization, etc. have been given by Pfeuty [10]. 

Our aim is to confirm the scaling behaviour of the structure factors for this model. 
In the disordered phase, the Hamiltonian for the model can be written as 

^ = E(l - ^t) - A E (38) 

i <ij> 

where the cxf = 2S°' are Pauli operators and the second sum is over nearest neighbour 
pairs. The critical point [WJ lies at A = 1, and the 1-particle energy is 

E{k) = 2A{k), (39) 

where 

A(A;) = [l + A2-2Acos(A;)]^/^ (40) 
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so that the 'critical wavevector' is /cq = and the energy gap is 

A = 2(1 -A). (41) 

The 1-particle exclusive structure factors have been discussed by Hamer et al. 
Multiparticle expansions for correlation functions for the quantum XY model in one 
space dimension have been obtained by Vaidya and Tracy [llj. The transverse Ising 
model is merely a special case of the model considered by them (Section 2.2 of Ref. 
[TT] for t = 0, 7 — > 1, and h = 1/A). Hence one can obtain exact expressions for the 
1-particle contributions to the correlation functions 

C'^'^in) = {SSS:)o (42) 



as: 



Cf^(n) = (1 - A')^/'-^ r dk cos{kn)A{k) (43) 
on Jo 



Hence one finds 
Si^ik) 



A2)l/4 



AA{k) 

Sf,{k) = \{l-X'y/'A{k) (44) 

In the vicinity of A ^ 1, A; ^ 0, equation fl39|) reduces to 

E{k) ^ A/(cg/A) (45) 

where 

fix) = /{! + x') (46) 



with c = 2, which is the expected relativistic form. The 1-particle structure factor Sfp 
reduces to 

5-(fc) ^ ^A-3/4|.'-(cg/A) (47) 

which has the expected scaling form (c.f. equation fl36l) ). with d = 1,-2 = 1,?7 = 1/4, i/ = 
1, the transverse Ising model values, and 

^[f{cq/A) = l//(cg/A). (48) 

The other transverse structure factor 

Sf^ik) - U'/%yyicq/A). (49) 

where 

l>;^/(cg/A) = /(cg/A). (50) 

Note that whereas Sf^{k) diverges as {A ^ 1, A; = 0}, Sfp^k) does not, and has a sub- 
leading critical index, two powers of A smaller than S^^. It appears that S^^ decouples 
from the one-particle state at the critical point. 
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The quasiparticle residue for the dominant spectral weight S^^ at A; = is 

Aik) = (1 - X^y/' [2(1 - A^l, (51) 

in agreement with Sachdev's result [1], after one takes into account differing 
normalization factors in our definitions. Note that in this case A{k) is independent 
of fc. 

We may deduce the scaling form of the full 1-particle structure function in the 
vicinity of the critical point: 

^-7/4 



-6{uj/A - /(cg/A)) 



(52) 



2/(cg/A) 

whence the scaling function for the dominant component may be taken as 

^^,;{cq/A,ujA) = ^r;icq/A)6{uj/A - /(cg/A)) (53) 

with normalization factor Zf^ =1/2. These are the simplest possible free-particle forms, 
save only the renormalization of the residue function with coupling. 



4. Comparison with Numerical Results 

4.I. The Transverse Ising model in higher dimensions 
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Figure 1. (Color online) 

Sfp (k) along high-symmetry cuts 
through the Brillouin zone for 
the transverse Ising model with 
coupling A = 0.1, 0.2, 0.3, 0.328, 
0.3284 on the square lattice. 
(From ref. [S]). 



Figure 2. (Color online) S'fp (k) along high- 
symmetry cuts through the Brillouin zone 
for the transverse Ising model with coupling 
A = 0.1, 0.2, 0.3, 0.328, 0.3284 on the square 
lattice. (From ref. [5J). 
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The behaviour of the transverse Ising model in higher dimensions is quahtatively 
similar. The 1-particle structure factors for the transverse Ising model on the triangular, 
square, and cubic lattices have also been calculated by Hamer et al. [5], using high-order 
series expansions. Some sample results for the square and cubic lattices are shown in 
Figures [TH 

For the square lattice, the critical point is estimated [12] to lie at A = 0.32841(2), 
and the critical exponents are expected to be the same as those of the classical 3D 
Ising model, namely rj = 0.0364(5), u = 0.6301(4), from various estimates [13]. The 
results for S^^ and Sfp along high-symmetry cuts through the Brillouin zone for the 
system with couplings A = 0.1, 0.2, 0.3, 0.328 and 0.3284 are given in Figures [H 
and [2j The results of a standard Dlog Fade analysis [5j of the series for (k) at 
k = (0,0) and k = {n/2,n/2) at k = (0,0), where the energy gap vanishes, give 
estimates Ac = 0.3284(4) with exponent —0.605(5), compared to the expected exponent 
iy{r) — 1) = —0.607. At momentum k = (7r/2,7r/2), where the energy gap remains 
finite, we find Ac = 0.34(3) with exponent 0.04(2) compared to the expected value 
urj = +0.0229. For Sf^, the estimate for the critical index is very close to the value 
z/(r/ + 1) = 0.65. 

In Figures [1] and [2] for A = 0.328 and 0.3284, we have biased the critical point to 
Ac = 0.32841 with critical index urj = +0.0229 in our analysis. We can see from these 
figures, that even for A = 0.3284 which is very close to the critical point, and Sf^ 
are still far from zero. This reflects the tiny value of the exponent rju, which implies a 
precipitous drop to zero just before the critical point. 

Figures [3] and H] show similar graphs for the simple cubic lattice. In this case, the 
critical point has been obtained previously [H] as Ac = 0.19406(6), and the critical 
exponents are expected to lie in the universality class of the 4D classical Ising model, 
where we expect the mean field exponents rj = , u = 1/2, modulo logarithmic 
corrections [1]. 

The analysis of S'J^p(k) at k = (0,0,0), where the energy gap vanishes, gives 
Ac = 0.19406(8) with exponent -0.54(1), while for Sf^{k) at k = (0,0,0), the estimate 
of the critical point is Ac = 0.194(4) with exponent 0.55(3). Away from k = (0,0,0), 
where the energy gap remains finite, we find Ac = 0.22(3) with exponent 0.03(2) for 
both 5'fp(k) and S'j'p(k). Allowing for logarithmic corrections, these estimates agree 
reasonably well with the expected values. 

In all cases, we see that the dominant structure factor Sf^ at the critical wavevector 
diverges at the critical coupling with exponent — z/(l — rj), while Sip vanishes with 
exponent consistent with z/(l + 77). Away from the critical wavevector, the structure 
factors both vanish at the critical coupling with a small exponent consistent with ut]. 
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Figure 3. (Color online) 
Sfp (k) along high-symmetry cuts 
through the Brillouin zone for the 
transverse Ising model with cou- 
pling A = 0.1, 0.15 and 0.19 on 
the simple cubic lattice. (From 
ref. 0). 



Figure 4. (Color online) Sfp(k) along high- 
symmetry cuts through the Brillouin zone 
for the transverse Ising model with coupling 
A — 0.1, 0.5 and 0.19 on the simple cubic 
lattice. (From ref. [5]). 



4-2. The Alternating Heisenherg Chain 

Schmidt and Uhrig [6] and Hamer et al. ^ have investigated the spectral weights of 
the alternating Heisenberg chain, which can be described by the following Hamiltonian 

H = Y. (S2, ■ S2.+1 + AS2.-1 ■ S2.) (54) 

i 

where the Sj are spin-^ operators at site i, and A is the alternating coupling. Here 
we assume that the distance between neighboring spins are all equal and the distance 
between two successive dimers is d. 

There is a considerable literature on this model, which has been reviewed by Barnes 
et al. [15]. At A = 0, the system consists of a chain of decoupled dimers, and in the 
ground state each dimer is in a singlet state. Excited states are made up from the three 
triplet excited states on each dimer, with a finite energy gap between the singlet ground 
state and the triplet excited states. This scenario is believed [T6l[T7l[18] to hold right up 
to the uniform limit A = 1, which corresponds to a critical point. At A = 1, we regain 
the uniform Heisenberg chain, which is gapless. 

Several theoretical papers |T9l[20l[2Tll22] have discussed the approach to the uniform 
limit. Analytic studies of the critical behaviour near A = 1 [19j have related the 
alternating chain to the 4-state Potts model, and indicate that the ground-state energy 
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per site eo(A), and the energy gap A (A) should behave as 

60(A) -eo(l)~5^/VlMV'5o)| (55) 
A(A) ~ 6'/y,/\Hm>}\ (56) 

as A ^ 1, where 6 = (1 — A)/(l + A). This corresponds to critical exponents a = 2/3, 
u = 2/3. The logarithmic terms in (135]) are due to the existence of a marginal variable 
in the model. 

For the uniform chain A = 1, and near kd 27r, Affleck |23| has obtained 
expressions for the correlation functions in the model, including logarithmic corrections, 
which correspond to an exponent rj = 1: 

1 (In r) 

G'ir) = mr) (57) 

Fourier transforming, one obtains the asymptotic form for S~~^{kd) as 

S{kd) ^ S-^{kd) = ^^^^\ ln(7r - kd/2)\'/' (58) 

Note that in this case (1 —ri)v = 0, so there is no power-law divergence in the structure 
factor, but rather a logarithmic one. 

This implies that for kd = 27r and as A — > 1, the asymptotic form for S{27t) diverges 

as 

S{27t) oc [- ln(l - X)f^ A ^ 1 (59) 

For < kd < 271, one expects S to be finite for any A. 

The results obtained by Hamer et al. [7J for S versus momentum k for A = 0, 0.6, 
and 1 are shown in Fig. O Note that S{k)dk = 2tt (here we set d = 1), independent 
of A, so the area under each curve is the same. Also shown in the figure are the results 
for S' = GSf— 27rln(l — ^)/fc]~^/^ at A = 1. The results appear reasonably consistent 
with the expected behaviour. 

For fixed values of k, Fig. [6] shows the integrated structure factor S versus A, where 
for each value of k, about 20 different integrated differential approximants to the series 
are shown. We can see that the results converge very well out to A = 1. The logarithmic 
divergence as A ^ 1 for the case kd = 2tt is clearly evident. 

For Q < kd < 27T, an analysis of the series for the 1-particle structure factor 
^ip = 5*1^^ using Dlog Fade approximants by Schmidt and Uhrig [6] appeared to show 
that it vanishes with a behavior close to (1 — A)^/^. Since S remains finite, one would 
thus expect that Sip/S vanishes like (1 — A)^/^. This agrees with a heuristic argument 
[6j that the 1-particle spectral weight should vanish like \/A, i.e. like S^^^/\ \n{5 /6o)\^^^, 
where 5 = (1 — A) /(I + A). It disagrees, however, with what one might expect from 
the transverse Ising model example, that the one-particle residue should vanish with 
exponent rjv = 2/3 at all wavevectors, leading to a behaviour (1— A)^/^. It is possible that 
a logarithmic correction term may again be disguising the true power-law behaviour; or 
alternatively, the power-law behaviour of the renormalized one-particle residue function 
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Figure 5. The integrated 
structure factor S versus k for 
the alternating Heisenberg chain 
with A = (dotted hue), 0.6 
(dashed hne), 1 (crosses). Also 
shown is the quantity S' = 
65'[-27rln(l-^)/A;]-3/2 f^j, ^ ^ 

1 (squares). (From ref. [7]). 



Figure 6. The integrated structure factor S 
versus A for the alternating Heisenberg chain 
with kd — 7r/2, tt, 3tt/2 and 27r. (From ref. 



might indeed be different away from the critical wavevector. It would be useful to have 
some further analytical guidance in this case. 

Fig. [7| shows numerical values from Hamer et al. [7] for the relative 1-particle 
weight Sip/S versus A at selected values of kd. It can be seen that for any non-zero 
value of k, Sip/ S decreases abruptly to zero as A 1. Only at kd = 0+, does Sip/ S 
remain finite (about 0.993) in the limit A = 1; but by then S has itself decreased to 
zero. 

Finally, we discuss the results for the spin auto-correlation functions, defined as 

1 poo 

$(cu) = — / dkS-+{k,uj) (60) 

Schmidt and Uhrig [6] argued that the critical behaviour for the total auto 
correlation function (summed over uj) of the 1-particle state $ip should be 



$ip oc (1 - A) 



1/3 



(61) 



modulo logarithms, as for the structure factors. 




Figure 7. The relative 1- 
particle weight Sip/S versus A for 
the alternating Heisenberg chain 
with kd = 0+, 7r/2, tt, 37r/2 and 
27r. (Fromref. |7]). 



Figure 8. The auto correlation functions 
versus A of the bilayer Heisenberg model for 
the 1-particle state ($ip), 2-particle states 
($2p), and two particle bound states Ti and 
T2. (Fromref. [8J). 



Figure [8] shows various auto-correlation functions versus A, reproduced from Hamer 
et ai. One can see that •I'lp vanishes at the hmit A = 1, while (1 — A)~^/^$ip increases 
almost linearly as A increases. The curve for ($ip + $2p); if we assume it is non-singular 
at A = 1 (i.e. the singularities in $ip and $2p cancel exactly), runs almost fiat with 
A once we neglect unphysical and defective approximants: that would indicate that 
the 2-particle sector accounts for about 99.8% of the weight, even at A = 1, which 
agrees almost exactly with the conclusions of Schmidt and Uhrig [6|. Remarkably, 
this is much higher than the fraction of 73% for the two-spinon continuum at A = 1 
calculated by Karbach et al. |24] from the exact solution. Also shown in Fig. [8] is the 
direct extrapolation of the 2-particle auto-correlation $2p using integrated differential 
approximants. These extrapolations assume that there is no singularity in $2p at A = 1, 
and the results give a somewhat smaller value of about 0.9 at A = 1. 

Overall, then, the 1-particle energy gap and spectral weight at general momenta 
appear to vanish as A — 1, following the behaviour predicted by Cross and Fisher [T9] . 
and already confirmed numerically by Singh and Zheng [25j. However, the 2-triplet 
spectral weight remains finite in the uniform limit and, in fact, appears to form the 
major part of the total spectral weight. Schmidt and Uhrig [6] already pointed out that 
indeed the 2-triplet states carry a larger portion of the total spectral weight than the 
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2-spinon states, calculated by Karbacli et al. [.24J. This argues that a description in 
terms of triplons remains equally valid with a description in terms of spinous for the 
uniform chain. 

4.3. Heisenberg Bilayer Model 

As our final example, we consider the Heisenberg bilayer antiferromagnet on the square 
lattice, with Hamiltonian 

H=JiJ2 E Sii-Sij + J2$:Sii-S2i (62) 

I = 1,2 <i,j> i 

where I = 1,2 labels the two planes of the bilayer. The physics of the system then 
depends on the coupling ratio A = Ji/J2- At A = 0, the ground state consists simply 
of 5* = dimers on each bond between the two layers, and excitations are composed of 
S = 1 'triplon' states [6j on one or more bonds. At large A, where the Ji interaction 
is dominant, the ground state will be a standard Neel state, with S = 1 'magnon' 
excitations. At some intermediate critical value Ac, a phase transition will occur between 
these two phases. It is believed that this transition is of second order, and is accompanied 
by a Bose-Einstein condensation of triplons/magnons in the ground state. 

Figures [9] and [10] show some series results for structure factors in the dimerized 
phase, calculated by Collins and Hamer [8]. Figure [9] shows the total static transverse 
structure factor ^(k) = S'"'"~(k) as a function of k at various couphngs A = Ji/J2- All 
results are for kz = vr, probing intermediate states antisymmetric between the planes, 
and we only refer to k = (A;^., ky) hereafter. 

The dominant feature is a large peak at the Neel point k = (7r,7r), which appears 
to become divergent as A ^ Ac, as we would expect. Figure [TO] shows the ratio of the 
1-particle structure factor S'ip(k) to the total ^(k) as a function of k. The 1-particle 
contribution generally remains the dominant part of the total, particularly near the Neel 
point. 

Let us now compare these results with theoretical expectations. From scaling theory 
(Sec. [2]), both the 1-particle structure factor and the total structure factor in the vicinity 
of the critical point should scale like (Ac — A)*^''"^-''^, at the critical (Neel) momentum. 
We expect this transition to belong to the universality class of the 0(3) model in 3 
dimensions, which has critical exponents [26] u = 0.707(4), 77 = 0.036(3), hence we 
expect {rj — l)v = —0.682(5), which is quite compatible with the numerical estimates. 

How does Sip behave at the critical coupling away from the Neel momentum? Here 
the behaviour is quite different from the previous models. The ratio Sip/ S decreases 
smoothly towards the critical coupling, and shows no sign of vanishing there. In fact the 
1-particle structure factor remains dominant everywhere, remaing at 80% of the total 
or more. Thus it appears that in this case the renormalized residue function does not 
vanish at Ac, except at the Neel momentum. 
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Figure 9. The total static 
structure factor >5'(k;) in the 
bilayer Heisenberg model as a 
function of k at various couplings 
A = Ji/ J2. (From ref. 0). 
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Figure 10. The ratio 5ip(k)/5(k) of the 1- 
particle static structure factor to the total 
static structure factor as a function of k 
in the bilayer Heisenberg model, for various 
coupHngs A — J\l J2- (From ref. [5]). 



5. Summary and Conclusions 

This paper consists largely of a review of the behaviour of structure factors near a 
quantum phase transition, at temperature T = 0. We have focused here on quantum 
spin models, but the conclusions should apply more generally. 

Section [2] reviewed current theory on the subject, drawn largely from Sachdev [1]. 
The generic scaling behaviour of both the total structure factor and the 1-particle 
exclusive structure factor is predicted to be the same, determined by the critical 
exponents r] and v. 

We then reviewed calculations of the structure factors for some specific quantum 
spin models. For the transverse Ising model in one dimension, exact results can be 
obtained [5] ; while for the transverse Ising model in higher dimensions [5] , the alternating 
Heisenberg chain [6l [7], and the bilayer Heisenberg model p], we have used some 
numerical results obtained from series expansions to high orders. For the most part, the 
results conform to theoretical expectations. 

Some significant differences have been noted, however, in the detailed behaviour of 
these models, particularly as regards the 1-particle structure factor. In the transverse 
Ising model the 1-particle residue vanishes at the critical point for all wavevectors, 
and so the 1-particle contribution to the total structure factor becomes negligible. For 
the solvable case of the one-dimensional chain, the residue is actually independent of 
wavevector. 

For the alternating chain, the one-particle residue again vanishes at the critical 
point, and it is the 2-particle 'triplon' state which appears to become dominant at the 
phase transition [6l [7]. But the residue appears to vanish with a different exponent 
depending on the wavevector, namely 2/3 at the critical wavevector and 1/3 away from 
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it, which seems pecuhar. It could be that the true exponent is disguised by logarithmic 
corrections, or perhaps the renormalized residue function does indeed behave differently 
at different wavevectors, and vanishes with a subdominant exponent away from the 
critical wavevector. Further analysis is needed here. 

For the bilayer Heisenberg model, on the other hand, the renormalized 1-particle 
residue vanishes at the critical wavevector only, and the 1-particle state remains 
dominant at the critical point. This is presumably the more typical pattern of behaviour. 
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